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From:  Saul  Blumenthal,  Ohio  State  University,  OSU  Research  Foundation  Project 

No.  716366 

RESEARCH  OF  SAUL  BLUMENTHAL  ON  AFOSR  GRANT  87-0306 

I.  ESTIMATION  WITH  TRUNCATED  DATA 

1.  Estimation  of  Sample  Size  with  Grouped  Data,  submitted  to  J.  Statist.  Planning  & 
Inference  (with  R.  C.  Dahiya). 

2.  Estimating  N  with  Percentile  Grouped,  and  Truncated  Data  from  a  Scale  Parameter 
Distribution,  Communications  in  Statist.  Theor.  Meth.17.  1988,  pp.3607-3628. 

3.  Estimating  the  Size  of  the  Sample,  to  appear,  Proc.  Indo-U.S.  Workshop  on  Bavestan 
Inference  in  Statist,  and  Econometrics.  1989,  Springer  Verlag. 

We  are  given  N  random  variables  Xj,...,  with  common  density  f(x;0 ) 
where  both  N  and  0  are  unknown.  The  values  of  the  X’s  can  be  observed  only  if  they  lie 
in  the  interval  ( 0,  T).  The  problem  addressed  is  estimation  of  N  and  0.  Maximum 
likelihood  estimators  and  variants  of  these,  namely  Bayes  modal  estimators  are  derived  and 
asymptotic  properties  are  studied. 

In  [1],  an  additional  restriction  is  imposed,  namely  that  the  interval  ( 0,  T )  is 
subdivided  and  only  the  number  of  X's  falling  in  each  subinterval  is  reported,  with  all  other 
details  about  actual  values  being  lost.  General  techniques  to  analyze  data  in  this  form  are 
developed,  and.restrictions  on  the  cstimadon  procedures  are  derived  to  assure  the  stochastic 
consistency  of  the  estimator  of  N.  A  widely  used  iterative  computing  scheme  is  shown  to 
fail  the  consistency  conditions  and  several  alternative  schemes  which  are  consistent  are 
examined.  An  analysis  of  a  software  reliability  data  set  is  carried  out,  and  Monte  Carlo 
comparison  of  several  estimators  is  performed,  leading  to  one  of  them  being  recommended 
for  use. 

In  [2],  instead  of  grouping  at  fixed  time  points  and  observability  on  an  interval 
of  length  T,  it  is  assumed  that  only  a  fixed  subset  of  the  quantiles  of  the  X’s  can  be 
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observed,  the  largest  one  being  the  (n/N)th,  with  n  <  N.  Problems  of  the  type  considered 
in  [1]  are  examined  for  this  alternative  model. 

The  paper  [3]  is  a  review  of  recent  work  on  all  aspects  of  estimation  with 
truncated  data,  but  concentrates  on  applications  of  this  work  to  estimating  the  number  of 
undetected  errors  remain  in  computer  software  after  an  initial  testing  period.  Both  the 
papers  in  the  software  reliability  literature  which  lead  to  truncated  data  models,  and  the 
papers  discussing  the  statistical  analysis  of  these  models  are  reviewed. 


IL-BURN-ItL 

1.  Bum-in  with  Mixed  Populations.  Ohio  State  University,  Statistics  Department 
Technical  Report  No.  372,  August  1987  (with  U.Q.  Pan). 

2. Bum-in  using  a  sequential  screening  procedure.  Technical  Report  1987.14,  Dept  of 

Mathematics  and  Statistics,  Wright  State  Univ.,  December  1987  (  with  U.Q.  Pan). 

3.  A  Reliability  Model  for  Bum-in.  Technical  Report  1987.12,  Dept  of  Mathematics  and 
Statistics,  Wright  State  Univ.,  November  1987  ( with  U.Q.  Pan). 

Closely  related  to  the  problem  described  in  Part  I  is  the  bum-in  problem.  In 
its  simplest  form,  a  batch  of  M  items  contains  N  (unknown)  defectives  with  known  failure 
distribution.  A  sequential  stopping  rule  is  desired  for  removing  enough  of  the  defectives 
to  assure  that  with  a  specified  high  probability  the  reliability  of  the  remaining  items  in  the 
lot  after  bum-in  exceeds  some  given  lower  bound.  Various  stopping  rules  have  been 
studied  for  both  exact  and  approximate  (large  lot  size)  properties.  The  technical  report  [1] 
is  based  on  the  dissertation  of  Dr.  Pan.  In  [2],  a  procedure  for  bum-in  based  on  a 
sequential  estimation  scheme  published  previorcW  by  the  proposer  is  extracted  from  [1], 
additional  computations  of  its  small  lot  size  propc  ni-  are  made  and  some  additional 
analytic  properties  are  developed.  This  report  is  cunently  being  readied  for  submission 
for  possible  publication.  In  [3],  the  characterization  of  the  fixed  bum-in  time  procedure 
which  achieves  the  after  bum-in  reliability  goal  is  extracted  from  [1].  This  fixed  time 
deper  Is  on  knowing  more  than  is  generally  available,  but  serves  as  a  benchmark  to  which 
to  compare  the  sequential  schemes  to  see  how  much  extra  time  they  take  because  the  extra 
information  is  not  available.  This  too,  will  soon  be  ready  for  submission  to  a  journal. 
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^RELIABILITY  DEMONSTRATION 

1.  Reliability  Demonstration  for  New  Series  Systems,  in  preparation. 

2.  A  New  Approximation  for  the  Reliability  of  New  Series  Systems,  in  preparation. 

3.  A  Two-Stage  Procedure  for  Reliability  Demonstration  of  New  Series  Systems,  in 
preparation. 

Series  systems  exhibit  wear  out  if  the  components  tend  to  wear  out  as  they  age. 
Acceptance  sampling  criteria  are  usually  stated  in  terms  of  Mean  Time  Between  Failures 
(MTBF)  for  an  aged  ( or  equilibrium)  system.  However,  the  acceptance  test  is  often 
administered  to  a  newly  produced  system.  If  the  test  does  not  take  into  account  the  fact  that 
the  new  system  has  a  much  larger  MTBF  than  it  will  have  after  aging,  too  high  a  proportion 
of  poor  quality  systems  will  be  accepted.  In  [1],  we  have  developed  tests  which  take  this 
aging  phenomenon  into  account  They  are  based  on  a  Poisson  approximation  for  the 
distribution  of  the  number  of  failures  a  new  system  will  have  in  a  fixed  time  period,  and  on 
the  assumption  that  the  component  failure  distribution  is  known  except  for  the  value  of  a 
scale  parameter.  We  have  also  looked  into  the  robustness  of  the  O.C.  curves  of  these 
procedures  in  case  the  failure  distribution  is  a  gamma  or  Weibull  with  unknown  shape 
parameter. 

In  [2], we  have  studied  alternatives  to  the  Poisson  approximation  and  found  that 
for  moderate  size  systems,  a  binomial  approximation  fits  the  true  distribution  mere  closely 
than  the  Poisson  or  any  of  several  other  possible  approximations.  We  are  now  looking  at 
test  design  using  the  binomial. 

In  [3],  we  are  studying  a  two  stage  testing  procedure  for  gamma  failure 
densities  which  develops  an  estimate  of  the  shape  parameter  based  on  the  first  stage  data, 
and  then  uses  that  estimate  to  develop  an  appropriate  test  time  and  acceptance  criterion  for 
the  second  stage  so  that  a  desired  O.C.  can  be  achieved. 


